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Abstract. Using a recently introduced linear reformulation of the Kuramoto model 
of self-synchronizing oscillator systems [1], we study a new class of analytically 
solvable oscillator systems defined by a particular coupling scheme. We show that 
these systems have a logarithmic scaling law in the vicinity of the critical point, 
which may be seen as anomalous with respect to the usual power-law behavior 
exhibited by the standard Kuramoto model. 

1 Introduction 

Throughout the physical world there are examples of systems of coupled oscillators — Josephson 
junction arrays, swarms of flashing fireflies, cardiac pacemaker cells in a heart, etc. — that 
spontaneously synchronize under certain circumstances. The Kuramoto model stands out as 
one of the few analytically soluble models that can describe the spontaneous synchronization 
phase transition of these systems j2|3| . For these reasons, it has become the paradigm in the 
field. In addition, one may think of such self-synchronization as a type of nonequilibrium phase 
transition, which highlights the potential usefulness of the Kuramoto model as a way to study 
this latter phenomenon. However, using the traditional approach, analytic solutions are largely 
restricted to systems where the number of oscillators approaches infinity (the continuum limit) 
and the oscillator populations are uniformly and globally coupled. 

A recent paper introduced a linear version of this Kuramoto model [IJ and discussed the 
properties of the reformulation, in particular how it permitted solution of systems containing a 
finite number of oscillators. In this paper, we work throughout in the continuum limit, where 
the linear reformulation of the Kuramoto model makes it possible to solve exactly a new set 
of coupled-oscillator systems. We shall focus on these systems' interesting anomalous scaling 
properties about the phase transition. However, before deriving and discussing these scaling 
properties, we will lay the foundation. We will begin by demonstrating the mapping of the linear 
reformulation onto the Kuramoto model. Next, we will derive the synchronization condition of 
the linear version. We will then solve the synchronization order parameter and critical point 
exactly for a class of systems with a specific coupling scheme, which cannot be solved using 
the traditional approach to the Kuramoto model. Finally we will outline the derivation of the 
anomalous scaling properties of the linear reformulation and discuss possible causes of this 
unusual behavior. 

The linear model can be expressed as 

ip{LU,t)^{ioj-'^)il;{uj,t)+ f2{uj,uj')g{Lu')il;{uj',t)duj' (1) 



where f2{uj, uj') describes the coupling between pairs of oscillators with characteristic frequencies 
UJ and Lo' respectively and the dot denotes a partial time derivative; the phases of the complex 
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variable ^{uj,t) correspond to those of the system's oscillators, the synchronization properties 
of which we will investigate; and g{uj') is the distribution of their characteristic frequencies. 7 is 
a parameter fixed according to the system parameters such that the amplitude of ^{uj,t) goes 
to a steady state in the long-time limit. This allows the linear model to be mapped onto the 
original Kuramoto model in the synchronized region, since with the nonlinear transformation 
= i?(w)e'^('^^ we can write the real and imaginary parts of eq. ([T]) as 

/>oo 

R{u) = -jR{uj) + / Sl{uj, uo')g{u')R{u') cos[6'(tj') - e{uo)]duj' (2) 



and, if R{uj) goes to a steady state in the long-time limit, eq. ^ is simply the Kuramoto model 
with a generalized coupling 

Kiu;,cj') = n{u;,cj')^. (4) 

H(Ll>) 

We discuss below how 7 is set such that this mapping can occur. 

Having reformulated the Kuramoto model in terms of linear dynamics, we can proceed 
to apply to it the tools and techniques developed for solving linear problems. Indeed, the 
synchronization problem can be discussed in terms of the spectrum of the linear operator on 
the RHS of eq. ([1]). More precisely, let IC{uj, oj') = n{uj, uj')g{uj') — iujS{uj — uj') and assume that 
the Fredholm integral equation 



duj'K:{uj,uj')(j>„{uj') ^ fi^4)„{uj), creZ,E (5) 
has a mixed, discrete-continuum spectrum {/inj/ij}. Then a generic solution of ([T]) is given by 
V^c^,*) = ^an0„Me(''"-^)*-f /" 6(a)0,He(''"-^)*da, (6) 

with coefficients {a„}„gz, {b(a)}^^s_ determined by initial conditions. 

To determine 7, we first solve for the spectrum {^„,/i(j}. We then set 7 equal to the real 
part of the eigenvalue with the largest real part. If there is only one eigenvalue whose real part 
equals 7 then in the long-time limit goes to a steady state, the linear model maps onto 

the Kuramoto model, and there is full phase locking and synchronization (as defined below). 
Otherwise, R{uj) will never go to a steady-state value, and the phases of "0 will never converge. 



Comments The model ([3]) with coupling constants ^ can be seen as a quenched limit of the 
full stochastic Kuramoto model [4j. Though deterministic, our model is more general because 
it allows for arbitrary coupling constants, while the standard stochastic version is understood 
only for specific distributions of these constants, with strong monotonicity constraints. As we 
shall see, the freedom introduced by our model leads to interesting qualitative differences for 
the behavior of the system in the vicinity of the critical point. 



2 Analysis of the Kuramoto model under the linear mapping 

Unless otherwise specified, we shall hereafter assume that we are in the situation where R{u!) 
goes to a steady state, and thus eq. ([3]) becomes the Kuramoto model with time-independent 
coupling. 

We now define our use of the term "synchronization" . We restrict the analysis in this paper 
to systems with no partial population of drifting oscillators, i.e. the incoherent-to-partially 
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locked (usually referred to as the synchronization transition) and the partially locked-to-fuUy 
locked phase transitions occur at the same point [5] . We say our system is synchronized if the 
synchronization order parameter given by 



(7) 



goes to a nonzero steady-state value. As mentioned above, this happens at the critical point 
where the real part of more than one eigenvalue becomes equal to 7. So in the synchronized 
region where the real part of only one eigenvalue equals 7, 

b{L0) 



dujg{ijj) 



(8) 



where &(w) is the eigenfunction corresponding to that differentiated eigenvalue, Xn- 

Let us consider now a specific example of global coupling in the linear model, J7(w, oj') = f2 
[B]. The linear model describes this system as 

/oo 
g{Lu')'il;{uj',t)dLu', (9) 
-00 

which maps onto the original Kuramoto model in the steady state with the following coupling 
constant: 

Ki.,.')^n.^I^^^, (10) 
y (w' - ujry + 7"^ 

where tOr is the collective frequency of the synchronized state and is given by the imaginary 
part of Ajv, 3[AAr]. In the following analysis (unless otherwise specified), we will assume that 
collective frequency 3 [A at] is zero as it can be eliminated through a phase change in ip, which 
effectively puts our system in a rotating reference frame and does not affect the analysis. 

With this coupling scheme, we can solve the spectrum of the RHS of eq. ([9|) exactly. Fortu- 
itously, eq. in]) maps onto the problem of stability of the fundamental mode of the incoherent 
state in the Kuramoto model with global uniform coupling 7\ , and we can follow the same steps 
towards a solution. With 7 set as described above, the spectrum comprises a continuous line of 
eigenvalues A in the complex plane along 3?[A] = —7 (where 3?[A] denotes the real part of A) for 
any value of f2 and one eigenvalue Aat at the origin, which stands apart from the continuum of 
eigenvalues if J7 > He- As SI ^ i7+, 7 — > 0. By contrast, for SI < Sic the continuum of eigen- 
values lies along the imaginary axis and Xn becomes indistinguishable from the continuum, as 
shown in Figure 1; and since, in the steady state, the entire spectrum remains, it is clear that 
r = where SI < Sic- That (in the steady state of the system) r — > as i? — > i7+ indicates that 
at that point the system is in a splay state [S] , where the phases of the oscillators in the system 
are locked in a uniform distribution from to 2tt. Also, above Sic, as the coupling between 
oscillators increases, r grows continuously from r = 0, a behavior akin to a second-order phase 
transition. This will be shown below. 

One can derive the nontrivial behavior of r from eq. ([9]). Assuming tp{uj,t) ~ b{uj)e^'^^, eq. 
dH]) becomes 

Ajv6(w) = iujb{uj) + SIA (11) 
where A = g{ijo)h{u})dui , which is a constant that can be determined self-consistently. 
Treating A as a, constant, one can algebraically solve eq. (fTT|) for b{ijj) and then plug this 
solution for b(uj) back into the definition of A to arrive at the following formula that determines 

l = nr -Mi^du; (12) 

J -00 -^N - 

As before, we set 7 = 5R[AAr]. Assuming g{uj) is an even function and nowhere increasing for 
uj >Q, then eq. (fT2|) has at most one solution for 7 and 7 > [9j, which implies 
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Fig. 1. The spectrum of eigenvalues associated with the RHS of eq. Q when Q > [2c- The spectrum 
comprises a continuum of eigenvalues along —7 and a single eigenvalue at the origin. 
7 — !■ 0, and the continuum approaches the imaginary axis and the eigenvalue at the origin. When 
J7 < f2c, 7 = 0, the continuum of eigenvalues sits on the imaginary axis and the eigenvalue at the 
origin becomes indistinguishable from the continuum. 



By taking 7^0+ (r ^ as discussed above) and assuming g(Lo) has a finite width, it becomes 
clear from eq. (fT3|) that 

= (14) 

Knowing that b{uj) = ^ ^^A^ from eq. pT|) . for Q > Qc (i-e. r > 0) we can also compute an 
explicit expression for r from eq. ([8]) (recalling that gioj) is even) to arrive at 

dwg{i^)—j= ^ (15) 

where 7 can be determined from eq. p3)) . So, for a given distribution g{Lo), eqns ^3]) . p4)) . and 
(fT5| completely specify r{fl) and Qc for K{uo,uo'), eq. pO|l . 



2.1 Scaling properties of the order parameter and coupling constant near the critical point 

Equations (jl3p and (|15p give interesting, new behavior (resembling a second-order phase tran- 
sition) in the critical regime. In order to study eq. (jl3p perturbatively, we start from its form 
([T^ and the identity d\^j^i^f{\isi ~iuj) = 0, which is valid for any smooth function /. Denoting 
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(in a distribution sense, assuming that we work in the correct Sobolev space [10] ) 

J — lU! 

we can expand around the origin to obtain 

Since 'Do{uj) — ■jtS{lu) (c.f. (fTi)) ). using the identity mentioned above, (fT^ becomes 



1 = n 



ng{0) - 3? i7 



'{co) 



doj 



for 7 — > 0, 



/-oo 7 - «^ 

where (/'(w) is the derivative of g^cu) with respect to w. This can be expressed as 

ujg'{uj) 



[2{ 7r.g(0) + 7 hm 



7^ + cj2 



duj 



which finaUy yields 



l = n[ ng{0) + 7 



duj I . 



Perturbatively, the behavior of 7 as i? ^ J7c becomes 



715(0) 



f2- f2c 



O 



[2- f2c 



(16) 



(17) 



(18) 



(19) 



(20) 



(21) 



where we have assumed the distribution g{Lo) is such that dug'[uj)/Lo is nonzero and finite. 

Let us now turn to eq. (fT5| . Using the Weber formula for the Laplace transform of the 
Bessel function of first kind, 



+ 72 



e-'^Vo(7s)ds 



we obtain 



r = 7 / g{s)Jo{js)ds, 



(22) 



(23) 



where g{s) = m{—s) is the Laplace transform of g{x), or the moment-generating function of the 
negative argument —s. Alternatively, we may say that r is proportional to the Hankel transform 
of g{s)/s. Eq. ((23)) allows us to extract the behavior of r(7) as 7 — > 0. First, notice that, since 
we are interested in the limit of vanishing 7, we may consider the test p.d.f.'s with compact 
support ghi'jj) = U[0, L] from the uniform family. Then the Laplace transform becomes 



1 - e 



-Ls 



9l[s) 



Ls 



I'iy), y = Ls, 



(24) 



where we introduce the scaled variable y = Ls. Clearly, the scaled function u has the behavior 

v{y)-h for y>l, 

(25) 

v{v) 1, for J/ < 1, 

so that eq. (|23p becomes 



27 

^(7) = Y / '^iy)Jo{'^y/L)dy, 



(26) 
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where the factor of 2 arises because we take the uniform distribution from —L to L. Let be 
the smallest zero of Jq, Jo(a;*) = 0. The integral eq. (|26p may then be approximated as 

where e is an infrared cut-off reflecting the small-y behavior of v eq. (j25p . Asymptotically in 
the small parameter 7/L, this gives the expansion 

r(7)- -25(0)7 log[.g(0)7] + O[5(0)7], 7^0, (28) 

where 5(0) = It is also worth mentioning that r^-y) is a monotonicaUy increasing function 
taking values from r(0) = to r'(7 — > 00) — 1, so it is a properly defined cumulative distribution 
function on M+. 



2.2 Specific examples of distributions 

With these general solutions for the parameters of any system with global coupling f2(w, to') = 
f2, we can solve for a specific system given its characteristic frequency distribution. Take for 
instance the Lorentzian distribution about LOr, i.e. g{u! — Wr) — TT[A'^+(t-u )^] • From eq. (jl4p we 
find i7c = A, and from eq. (fT2)) . 7/0^ = i7 — Zi. Using these in eq. p^)) . we obtain 

2C0S-1'' 



Hor = , = (29) 



ttWI 



for J7 > i7c, as shown in Figure 2. ?■ = for ^? < i7c as discussed above. From a Taylor expansion 
around the origin of this expression, we obtain the following scaling for this distribution: 

This agrees with eqns (I?!]) and (|28p knowing that, for this Lorentzian distribution, g{0) = ^ 
and dujg'{uj)/ui = 

Similarly, for a uniform distribution about ujr, i.e. g(cij — 0;^) = for |a; — oj^l < TrA/2 and 
otherwise, the above equations give Qc = A, ^unif = ^ cot (f^), and 



2 ^"''V2fi/ 

tan 

2 i7 



(31) 



for fi > flc (see Figure 2). As above, r = for J7 < Qc- Using a Taylor expansion about the 
origin, the scaling for the uniform distribution is then 

Again, there is agreement with eqns (I^T]) and (|28p as (/(O) A- and dtog' {ijj) /uj = — (T^p-- 



3 Discussion 

The synchronization phase transition discussed in this paper is of the nonequilibrium variety, 
but there are many interesting connections with the more typical equilibrium phase transitions. 
One way to illuminate the nature of this nonequilibrium phase transition and create a platform 
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Fig. 2. The synchronization order parameter as a function of the normalized coupling constant for a 
uniform and a Lorentz distribution of characteristic frequencies. 



for comparison with standard theories of equilibrium phase transitions would be to determine 
the critical behavior of the fluctuations of the order parameter close to the critical point, which 
exhibit power-law scaling in the uniform Kuramoto model . 

Although the physical reason for the logarithmic scaling behavior of the synchronization 
order parameter close to the critical point in the linear model eq. ([T]) remains unclear, we wish 
to highlight a number of distinctive features of this phase transition that may be a cause of 
such behavior: 

(1) The specific coupling considered in this paper (eq. 10) is of a very different nature from 
the standard uniform coupling scheme associated with the Kuramoto model. It is therefore 
unsurprising that the scaling behavior of the linear reformulation also differs from the usual 
square-root scaling behavior of the normal Kuramoto model close to the critical point. This 
particular coupling is asymmetric and has the interesting property that the coupling be- 
tween two oscillators becomes stronger the closer one oscillator's characteristic frequency is 
to LOr and the further the other's characteristic frequency is from ujr (reminiscent of problems 
involving asymptotic freedom [H]). It should also be noted that close to the critical point, 
i.e. 7 ^ 0, if the characteristic frequency of one of the oscillators is uJr then the coupling 
tends to infinity in one direction and zero in the other. 

(2) The standard Kuramoto model describes two phase transitions — one from incoherence to 
partial locking and a second from partial locking to full locking — which in general occur at 
different points in parameter space [S]. In our paper we consider systems where these two 
phase transitions are coincident, i.e. where a bicritical point occurs. 

(3) This logarithmic scaling near a critical point is interestingly reminiscent of the behavior of 
the marginal eigenvalues in standard renormalization group theory [13j . 
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Although this paper only explored the critical properties of one particular type of coupling, 
the full range of coupling schemes K{uj,uj') (eq. ([4])) that can be analyzed with this method 
remains to be determined. 
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